Impact of weak localization in the time domain 
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We find a renormafized "time-dependent difltusion coefficient", D{t), for pulsed excitation of a 
nominally diffusive sample by solving the Bethe-Salpeter equation with recurrent scattering. We 
observe a crossover in dynamics in the transformation from a quasi- ID to a slab geometry imple- 
mented by varying the ratio of the radius, R, of the cylindrical sample with reflecting walls and the 
sample length, L. Immediately after the peak of the transmitted pulse, D{t) falls linearly with a 
nonuniversal slope that approaches an asymptotic value for R/L 3> 1. The value of D{t) extrapo- 
lated to t = depends only upon the dimensionless conductance g for R/L <^ 1 and only upon kio 
for R/L ^ 1, where k is the wave vector and £o is the bare mean free path. 
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Weak localization (WL) of electronic and classical 
waves arises from the interference of counterpropagating 
partial waves in closed loops. Its impact upon electronic 
conductance has been widely studied using steady- 
state methods such as magnetoresistance 0, and can 
be directly visualized as enhanced retroreflection of light 
from random samples Q. Partial waves with trajectories 
over a wide range of lengths contribute to these phe- 
nomena. It is, however, of great interest to investigate 
the variation of WL upon pathlength 0, IE IE S since 
the impact of WL presumably increases with pathlength. 
Though it has not been practical to isolate paths of spe- 
cific lengths in studies of electronic conductance, this 
can be accomplished for classical waves in time-resolved 
measurements of pulsed transmission through random 
media H H E 111 El 111. Time-resolved transmis- 
sion measurements have generally been consistent 
with diffusion theory, which predicts a simple exponen- 
tial decay of the average transmission following an initial 
rise, with the decay rate due to leakage from the sam- 
ple of I/Id = n^Do/{L + 2zo)^, where to is the diffu- 
sion time, Do is the diffusion coefficient, L is the sample 
thickness, and zo is the extrapolation length. Recently, 
however, Chabanov el al. |3l observed nonexponential 
decay of pulsed microwave transmission in a quasi- ID 
sample which they characterized by a "time-dependent 
diffusion coefficient", D{t). The decrease in the decay 
rate with time was compared to the leading correction 
term in a supermatrix model calculation by Mirlin 7] of 
the tail of the electron survival probability. In the case of 
the orthogonal ensemble, it can be expressed in term of 
the renormalization of a "time-dependent diffusion coeffi- 
cient", D{t)/Do — 1 — (2/7r^(7<£))t, where g is the dimen- 
sionless conductance. A linear decay of D{t) was found in 
microwave experiments; however, D(t) did not extrapo- 
late to the bare diffusion coefficient at t = and the 
scaling behavior of the slope of D{t) did not agree with 
the above expression. Nonexponential decay of pulsed 
transmission has also been reported recently in numeri- 



cal simulations in 2D 0| and in a self-consistent diffu- 
sion theory, which includes recurrent scattering For 
times much larger than the Heisenberg time, tn = gtu, 
it has been predicted that the decay of the electron sur- 
vival probability follows a log-normal behavior due to the 
presence of prelocalized states . These long-lived states 
might be associated with rare configurations of disorder 
in the medium 

In this Letter, we solve the Bethe-Salpeter equation 
with recurrent scattering included in a self-consistent 
manner that satisfies the Ward Identity to obtain the 
average time-dependent intensity transmitted through a 
random sample following pulsed excitation. The sam- 
ple is cylindrical with refiecting walls. By changing 
the ratio of the longitudinal dimension L and the ra- 
dius R, a smooth transition can be made between two 
key experimental geometries: a quasi-lD geometry with 
L ^ R, which is commonly employed in microwave ex- 
periments, and a slab with L <^ R, which is the typical 
optical geometry. For the quasi-lD geometry, we find 
D(t)/Do = A - {2B/TT'^gtD)t for t <€. tn- The constant 
part. A, is universal, depending only on g, while B is 
nonuniversal and depends upon L/Iq, where £o is the bare 
mean free path, as well as upon g. In the limit g ^ 1 
and L/£q 3> 1, our results coincide with supersymmetry 
calculations by Mirlin 0. For intermediate values of 
g > 5, our results are in agreement with experiment 
For the slab geometry, D approaches a nearly constant 
renormalized value, being equal to Do{l — 1.03/(/c£o)^) 
in the limit L /£o ^ 1 . This is close to the result of WL 
theory for a bulk system. 

We consider a scalar wave incident on the front sur- 
face of the random sample at z = 0. We assume that the 
medium possesses neither absorption nor gain and that 
scattering is isotropic. The time evolution of intensity 
within the sample is obtained from the Fourier trans- 
form of the frequency correlation function Co(a;,r) = 
{(j)Q+{r) (r)), where il* — nzLuj/2, fl is the wave fre- 
quency, Lu is the frequency shift, and (/'o(r) is the wave- 
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function or field inside the sample [lol . The function 
Co(w, r) is obtained from the space- frequency correlation 
function Co(a;,r,r') = (0^;+ (r) (/)^_ (r')), which satisfies 
the Bethe-Salpeter equation, 

Co(^,r,r') = (0,,.(r))(0*_(r')) 

+ J dn dT2 dr3 dri{Gn+ (r, Yi)){Gn- (r', rg)} 
X [/r!(t^;ri,r2;r3,r4)Cr2(tJ,r2,r4), (1) 

where ((/)n(r)) represents the coherent source inside the 
sample and (G'n(r,ri)) = — ^^^^jjzt^j^ is the ensemble- 
averaged Green's function. Here K — k + where k — 
^ , u is the wave speed, £ is the scattering mean free path 
determined from the imaginary part of the self-energy of 
(G) . £o is determined from the single-scattering diagram 
only. The vertex function Un in Eq. (1) represents the 
sum of all irreducible vertices. Here we approximate Un 
as 

t^a(w;ri,r2;r3,r4) = 

^ [1 + Siu;, k)] <5(ri - V2) S{r, - rg) 6{n - r4). (2) 

The first term, proportional to 47r/£o, generates all self- 
avoiding multiple-scattering diagrams, which produce the 
diffusion result when L £o 18], whereas the factor 
d{Lu,k) in the second term represents the WL contri- 
bution to Uq. The presence of the second term in the 
vertex function renormalizes the mean free path, giving 
^(w, fc) = fo/[l + ^(w, fc)]. The Ward Identity, which en- 
forces flux conservation, requires that the same £{lo^ k) 
should appear in (G) of Eq. (1). Here, we assume that 
the system is far from the localization threshold and that 
the renormalized mean free path is scale independent. In 
a bulk system, (5(0;, fc) can be obtained from the renor- 
malized frequency-dependent diffusion coefficient |l7l |. 

where D{u),k) — v£{lu, k)/3, and G is the Green's func- 
tion for the diffusion equation and refiects the return 
probability. Here, we solve for G by using the boundary 
conditions of a reflecting tube. After taking the spatial 
average, we find 



with 

and 
62(0;, k) 



S{uj, k) = 5i{uj, k) + (52(^, fc), 



V 1 
di(uj, k) = ^ > — : — - 

71—1 "' 



(4) 
(5) 



where N = [kRY /A is the number of transverse modes, 
qn = mT/L,^nc = aL/irlo, and L = L + 2zo with 
zq = 0.71^0 Um- I'^ the above equations, a/£o denotes 
the upper momentum cutoff, q,., and is chosen as l/£o 
in our calculations unless specified otherwise. The i5i 
term arises from diffusive modes which are uniform in the 
transverse directions, i.e., q± = 0. When L/£o is large, 
the static limit of this term becomes 6i{uj = 0) = 1/3(7, 
where g — AN£o/3L. This term is responsible for the 
linear decrease of the decay rate in the transmitted in- 
tensity {I{t)). The 62 term represents the contribution 
from other transverse modes below q±^c- This term be- 
comes important only when R ^ £q and is responsible 
for the renormalization of diffusion constant in the limit 
of a slab geometry. In our calculation, we replace Dq in 
Eqs. (5) and (6) by v£{uj,k)/3 and solve for S{uj,k) in a 
self-consistent manner. 

For simplicity of calculation, we assume both the exci- 
tation intensity and the scattered intensity are uniform 
in the transverse cross-section of the tube. Eq. (1) can 
now be written as 1131 



, lUJZ z 

Go(w,z) = exp — 

V £0 



1 + S{k,iu) 
Ati£q 



dz'Hiu;,z~ z')CniuJ,z'), (7) 



where 



(tj, z — z') = -K I dp' 



exp 



( V - + (z - z') 



(6) 



p2 + (z - z'Y 

Eq. (7) is solved numerically for Cq{uj,L). Its Fourier 
transform gives {I(t)), from which we can calculate 
the renormalized "time-dependent diffusion coefficient" , 
D{t)/Do = -to dln{I{t))/dt, for t > to ^. 

We first consider a quasi-lD geometry with L :s> R — 
£q. In this case, since 82 is nonzero only when R > TT£o/a, 
we can set 82 = 0. For the case of iV = 800, D{t)/Do 
is plotted in Fig. 1 as a function of the dimensionless 
time, t" = 2t/TT'^gtD, for L/£q=2Q, 40, and 80. Also plot- 
ted in Fig. 1 (dotted line) is Mirlin's analytical result 
■7], D{t")/Do = 1 - t". The decay rate of D{t")/Do 
appears to decrease with increasing LI£q and is always 
smaller than the analytical result. If we fit the lin- 
ear region to D{t")/Do ^ A — Bt" , the intercept A is 
found to decrease with increasing L/£q. These observa- 
tions are consistent with recent microwave experiments 
1^. In order to understand the behavior of A and B, 
we have carried out a systematic study by varying L/£q 
and N . The results for A are shown in Fig. 2, where 
A is plotted as a function of g for L/£q^1{), 20, and 
40. All the data points can be well fit by a single curve, 
A ^ 1.00 - 0.27/g - 0.17/g2, indicating that ^ is a func- 
tion only of g. In the limit of large we recover Mirlin's 
result, A = 1. It is interesting to note that the above re- 
lation gives, 13 = dlnA/dliiL = -(O.27/5 + 0.34/g2)/A, 



.(8) 
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FIG. 1: D{t")/Do plotted as a function of the dimensionless 
time, t" = 2t/TT^gtB, for quasi-lD samples with A'' = 800. 
The dotted line is the analytical result of Ref [7], 1 — t" . 
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FIG. 2: The parameter A obtained from the fit of D{t")/DQ 
with A—Bt" is plotted as a function of g for quasi- ID samples. 
The dashed curve is ^ = 1.00-0.27/g- 
ci versus L/lo\ the dashed curve is ci 



0.17/g . In the inset: 
= 1.03 - 4.42^o/i- 



which corresponds to the exponent of the local scaling 
relation, A cx L'^. If we use the value in the microwave 
experiments of 5 = 5.5, we find j3 = —0.06, which is 
in good agreement with the measured value of —0.05. 

The corresponding results for B are shown in Fig. 3. 
Unlike A, B depends not only on g but also on L/^o- 
For each L/(.q, we fit B with Bq + Bi/g + Bijg^ , where 
Bi = Bi{L/£Q). The values of Bo are found to be linear 
in {L/£a)^^, as shown in the inset. In the limit of large 
L/io, Bo approaches the Mirlin's result of i? = 1. Thus, 
our calculations indicate that Mirlin's result is valid only 
when both g and L/£o are sufficiently large. In Fig. 4, 
we compare our theory with microwave measurements 
for iV = 68 and L/£o = 6.4 (Sample A), 9.5 (Sample 
B), and 19.2 (Sample C), corresponding to £0 — 9.5 cm. 
The best fit is found for a = 0.5, and good agreement 
is found for Samples A and B. For Sample C, for which 



FIG. 3: The parameter B obtained from the fit of D{t")/Do 
with A — Bt" is plotted as a function of g for quasi- ID 
samples. The dashed curves are the fit of B with B = 
Bo + Bi/g + B2/g^, where Bi are constant. In the inset: 
Bo versus {L/£o)~^; the dashed line represents a linear fit. 
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FIG. 4: D{t")/Do plotted as a function of t' = t/to for three 
samples studied in Ref. [8], with A'^ = 68, L/£o — 6.4 (Sample 
A), 9.5 (Sample B), and 19.2 (Sample C), and Do = 41.5 
cm^/ns. The theoretical results for the corresponding samples 
are shown as solid curves. 



(7 = 3, however, our theory predicts a smaller decay rate, 
reflecting the important role of localization effects, which 
are not included here. 

We next consider the crossover from a quasi- ID to a 
slab geometry, as R/L increases. When 2R/L is suffi- 
ciently large, the Si term falls leading to a slope of D{t) 
with a magnitude which decreases as l/g oc 1/i?^. In 
the limit of large R/L, the function D{t) approaches a 
nearly constant value, D, which is determined by the 62 
term. Thus, we focus our discussion on A. In Fig. 5, we 
present A as a function of g for different values oi R/£o 
for i/^o = 20. These results suggest, A = 1 — co(L, R)/g- 
From the fit (dashed lines) to the data, we find that 
co{L,R) is proportional to R^ when 2R > L, suggest- 
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FIG. 5: The parameter A obtained from the fit of D{t")/DQ 
with A — Bt" is plotted as a function of g for L/la — 20 and 
different radii R/lo- The dashed curves are the fit of A with 
A — Ao + A\/ g + A2/ , where Ai are constant. In the inset: 
A versus 2R/L for Mo = 4 (A), 6 (B), 8 (C), and 12 (D). 



ing A = \ — ci{L) / iklof' . This can also be seen in the 
inset of Fig. 5, where A is replotted as a function of 2R/ L 
for different values of Mq. It is seen that A approaches 
a constant value when 2R > L for each value of Mq. In 
order to find the behavior of ci(L), we repeat the calcu- 
lation for L/£o = 10, 30 and 40. The results are plotted 
in the inset of Fig. 2. The fitting of these data suggests, 
ci = 1.03 — AA2£o/L (dashed line), which is the behav- 
ior of the renormalized diffusion coefficient Z? in a slab. 
Thus D decreases with increasing L, which is due to the 
presence of longer recurrent scattering paths in a thicker 
sample. This behavior is different from that found in thin 
samples in optical measurements when L/Iq < 6 |12|. It 
is interesting to note that when L/£q ^ 1, our result 
D/Dq ~ 1 — 1.03/(/c£o)^ is close to the static result of 
Eq. (3) in the bulk for a = 1, i.e., 1 - 3a/7r(Mo)2. 

In conclusion, we have solved the Bethe-Salpeter equa- 
tion with recurrent scattering to obtain WL in the time 
domain. Following peak of the transmitted pulse, the 
"time-dependent diffusion coefficient" D{t) falls nearly 
linearly. We find the extrapolated value of D{t) at t = 
and the slope of D{t) for samples with different scattering 
strengths, as the aspect ratio of the sample is changed, 
transforming the sample geometry from a quasi- ID to a 



slab. From this prospective, WL in steady-state mea- 
surements can be understood in terms of the increasing 
impact of WL associated with trajectories of increasing 
length. 
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